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1.  INTRODUCTION 


The  accuracy  of  the  Navy  navigation  satellite,  Tran¬ 
sit,  is  dependent  on  precise  orbit  determination  and  orbit 
prediction.  Improvements  in  orbit  prediction  have  been 
limited  by  uncertainties  in  orbit  disturbance  from  solar 
pressure  and  atmospheric  drag.  The  Disturbance  Compen¬ 
sation  System  (Discos)  was  designed  at  Stanford  University 
under  subcontract  to  the  Applied  Physics  Laboratory  of  The 
Johns  Hopkins  University.  It  was  flown  on  the  Triad  satel¬ 
lite,  and  it  has  opened  possibilities  for  improvements  in 
accuracy  and  operational  convenience  for  the  Transit  sys¬ 
tem  by  freeing  Triad's  orbit  of  disturbances  larger  than 
5  x  10“12  g4 

To  eliminate  the  effect  of  forces  on  the  Discos  sys¬ 
tem  because  of  the  various  satellite  components,  each  main 
electronics  package  (shown  in  Fig.  1)  had  its  mass  proper¬ 
ties  measured.  To  accomplish  this  a  bifilar  pendulum  was 
designed  that  could  rotate  a  body  into  six  distinct  positions 
with  respect  to  a  coordinate  system  that  had  its  origin  fixed 
at  the  body's  center  of  mass.  With  this  rotational  capability, 
six  independent  moment  of  inertia  measurements  were  made. 
The  data  were  then  reduced  with  the  aid  of  a  computer,  and 
the  complete  inertia  matrix  with  respect  to  the  package  cen¬ 
ter  of  mass  was  determined. 


Preceding  page  blank 
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2.  ANALYSIS 


Figure  2  shows  a  line  with  respect  to  a  mass  parti¬ 
cle  mj.  The  moment  of  inertia  of  the  particle  with  respect 
to  the  line  L  (X,  /i,  y)  is 

=  X^Ixx  +  M2Iyy  +  y2Izz  -  2MXlxy  -  2Xyixz  -  2MVlyz  .  (1) 

X,  fi,  and  Oi  are  the  direction  cosines  of  a  line  L  in  space 
that  passes  through  the  origin  of  a  fixed  Cartesian  system 
of  coordinates. 

In  order  to  solve  for  the  unknown  cross  products  in 
Eq.  (1),  it  is  necessary  that  the  body  shown  in  Fig.  2  be 
moved  to  six  different  angular  positions  and  that  the  moment 
of  inertia  is  measured  at  each  of  those  positions.  Equation 
(1)  was  derived  so  that  only  three  of  the  angular  positions 
can  be  mutually  perpendicular.  Figure  3  shows  three  of  the 
mutually  perpendicular  axes  x7,  y7,  and  z7.  Each  of  these 
axes  could  be  placed  on  the  axis  of  rotation.  Also,  the  three 
skew  lines  are  shown  in  Fig.  3;  their  angular  rotations  and 
line  numbers  are  indicated  below: 

Line  Rotation 

1  +5° 

2  -5°, +45° 

3  +45° 

To  use  the  moment  of  inertia  data  obtained  from  the 
six  different  measurements,  the  equations  are  set  up  as 
shown  in  Eq.  (2).  Equation  (1)  is  rewritten  w.here  the  prime 
refers  to  the  moment  of  inertia  of  the  axis  being  measured: 
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Fig.  2  LINE  WITH  RESPECT  TO  A  MASS  PARTICLE 
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fiXlxy  +  Xylxz  +  M^Iyz  =  i[X2Ixx  +  M2iyy  +  y2izz  -  IT  ]  .  (2) 


Now,  setting  up  a  matrix  of  the  left-hand  side  of  Eq.  (2)  in 
consideration  of  Fig.  4  where  Sa  =  sin  ot  and  Ca  =  cos  a 
(also  for  three  skew  lines)  we  obtain  the  following: 


“S^Sot  C0  Sa  Q0  SotjCa  SjSjSo^Ca" 

! 

1 

►— < 

l _ 

S^2Sa2C^Sa2  C^S“2Ca2  S^2S“2C“2 

1 

Ixz 

Q«3s«8c»3 

iyz 

In  setting  up  the  right-hand  side  in  a  matrix  form,  we  ob¬ 
tain: 


(C8  Sa  )2Ixx  (SJ3  S«  )2Iyy  Caizz  -  1^  " 
(C^2Sa2)2l^x  (S82SO!2)2I^y  CQ^Izz  -  £ 
(C83SO!3)2I^x  (Sj33S«3)2Iyy  C^Izz  -  I* 


(4) 


Using  matrix  notation  on  the  left-  and  right-hand  side, 

Ixy 

I 

[Di  -  ILil  (})  ,  (5) 


rcij] 


Ixz 

Iyz 


therefore  letting 


Ixy 

Ixz 

iyz 


=  Eij.  Solving  for  the  cross  products 
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(Ixxlyylzz  -  Izylxx  -  Izxlyy  -  Iyxlzz  +  21yzlzxlyx) 

2  2  2 

+  Xdxxlyy  +  Iyylzz  +  Ixxlzz  -  Izy  +  Izx  +  lyx)  (8) 

+  X^(Ixx  +  lyy  +  Izz)  -  X^  =  0  . 

The  three  roots  of  the  above  will  be  the  moments  of 
inertia  about  the  principal  axis  of  the  package. 


Since  for  rotation  about  a  principal  axis  the  angular 
velocity  vector  coincides  with  this  axis,  the  set  of  numbers 
oiy,  w'z,  which  satisfies  the  following  corresponding  to 
Ij  of  Eq.  (8),  consists  of  the  direction  numbers  for  axis  I^. 

(Ixx  -  I  )co  ^  -  Ixyoo  ^  -  Ixzu)  ^  =  0 
lx  ‘  y  z 

-IyxcoJ1^  +  (lyy  -  I  )co  -  IyzW  ^  =  0  (9) 

x  1  y  z 

-Izxo)  ^  -  Izyco  ^  +  (Izz  -  I  )co  ^  =  0  . 
x  y  1  z 

Solving  Eq.  (9)  will  result  in  the  ratios  :  WyV^  : 
and  the  axis  associated  with  1^  is  thereby  defined  relative 
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to  the  original  coordinate  system.  To  solve  the  ratios,  let 
one  direction  number  equal  1„0.  Repeat  this  procedure  for 
the  other  two  roots  of  Eq.  (8)  to  define  the  other  principal 
axes  with  respect  to  the  original  coordinate  system. 

By  using  the  above  analysis  and  by  knowing  the  mo¬ 
ment  of  inertia  of  a  package  fo:.-  three  perpendicular  axes 
as  well  as  three  skew  axes,  tho  inertia  matrix  can  be  com¬ 
pleted.  Once  this  inertia  matrix  is  completed,  the  princi¬ 
pal  moments  of  inertia  of  the  body  can  be  determined  as 
well  as  the  orientation  of  the  principal  axes  with  respect  to 
the  coordinate  axis  system  usee  for  the  original  measure¬ 
ments. 
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3.  DEVELOPMENT  OF  THE  MEASURING  TECHNIQUE 


The  design  goals  for  the  measurement  of  the  mass 
properties  were  as  follows: 

1.  Error  in  the  moment  of  inertia  +0.  5%  of  the  sum 
of  the  three  inertias  about  any  set  of  three  mutually  perpen¬ 
dicular  axes; 

2.  Mass  center  location  determined  within  +0.  001 
inch  of  mounting  surfaces  and  hole  patterns;  and 

3.  Mass  measurements  of  all  packages  to  within 
±0.  050  gram. 

Of  all  the  requirements,  the  third  was  obviously  the 
one  that  would  be  the  easiest  to  meet.  To  determine  the 
inertia  matrix  and  center  of  mass  it  was  apparent  that  ex¬ 
traordinary  precision  would  have  to  be  obtained. 

There  were  two  basic  proposals  put  forth  by  APL  to 
obtain  the  complete  inertia  matrix: 

1.  Obtain  a  large  A-frame  and  suspend  a  torsional 
pendulum  at  least  6  feet  long  from  it,  and 

2.  Develop  a  pendulum  similar  to  a  trifilar  pendu¬ 
lum. 


The  second  proposal  was  selected  for  the  following 
reasons: 


1.  Package  orientation  and  subsequent  dynamic  un¬ 
balance  would  have  a  minimal  effect  on  the  system. 

2.  A  preliminary  calculation  of  the  error  allowed 
on  the  smallest  package  showed  it  to  be  on  the  order  of 

2.  5  x  10"5  slug  ft2.  It  was  apparent,  therefore,  that  air 
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drag  on  any  system  used  would  affect  the  period  measure¬ 
ment  and  should  be  eliminated  if  possible.  By  using  the 
second  proposal  the  entire  system  could  be  put  efficiently 
into  a  bell  jar. 

3.  Other  considerations  would  be  the  weight  and, 
therefore,  the  configuration  of  the  basic  system.  A  bifilar 
pendulum  lends  itself  to  efficient  mass  distribution.  Most 
of  the  mass  can  be  concentrated  at  the  center  of  oscillation, 
thereby  eliminating  large  Mr2  terms. 

Figures  5  and  6  show  the  preliminary  bifilar  pendu¬ 
lum  arrangement.  Part  of  the  testing  on  this  bifilar  pendu¬ 
lum  concerned  the  problem  of  determining  the  best  suspen¬ 
sion  system  wire  diameter  and  material.  The  figures  show 
the  pendulum  supported  by  two  0.  007  inch  BeCu  wires.  Sub¬ 
sequent  to  the  use  of  the  0. 007  inch  BeCu  wire,  music  wire 
and  hypodermic  needle  tubing  suspension  systems  were 
tested.  The  final  suspension  system  chosen  was  stainless 
steel  hypodermic  needle  tubing  with  a  0. 042  inch  O.  D.  and 
a  0. 006  inch  wall. 


The  basic  equation  for  determining  the  moment  of 
inertia  on  a  bifilar  pendulum  (knowing  the  period)  is 


t2wd2 

le^L 


(10) 


where  I  =  moment  of  inertia,  T  =  period,  D  =  distance  be¬ 
tween  the  supporting  wires,  W  =  weight  on  the  pendulum, 
and  L  =  wire  length. 

Figure  6  shows  the  pendulum  with  a  test  mass  in 
place  and  small  pins  that  could  be  selectively  removed  and 
relocated  to  produce  a  change  in  its  moment  of  inertia. 
This  test  mass  proved  that  a  moment  of  inertia  change  of 
approximately  +0. 5%  of  the  sum  of  three  perpendicular 
axes  could  be  measured.  After  it  was  proven  that  the  bifi¬ 
lar  pendulum  would  give  the  accuracy  desired,  the  design 
of  a  more  sophisticated  system  was  begun. 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

Silvi*  Sp*'*6  Maryland 


4.  FINAL  PENDULUM  DESIGN 


Based  on  the  analysis  and  the  preliminary  experi¬ 
ments,  a  platform  was  designed  that  could  rotate  to  three 
perpendicular  axes  and  that  could  also  be  offset  to  locate 
the  package  into  three  skew  positions.  Figures  7  and  8 
show  a  package  in  two  of  the  orthogonal  axis  positions,  and 
Fig.  9  shows  a  package  in  one  of  the  skew  positions.  A 
counterweight  was  mounted  beneath  the  package  on  the  plat¬ 
form  frame  to  compensate  for  the  unbalanced  moment 
caused  by  the  location  of  the  package  in  its  various  posi¬ 
tions. 

Inherent  in  the  platform  design  was  the  necessity  to 
design  it  so  that  the  package,  when  attached  to  the  platform, 
could  be  positioned  at  its  center  of  mass.  In  Fig.  7  the 
arrows  show  the  adjustment  methods  as  indicated  below: 

Function  (see  Fig.  7) 

This  shaft  was  adjusted  by  a  screw  on 
the  bottom  that  allowed  the  package  to 
be  adjusted  in  the  vertical  direction. 

This  shaft  was  adjustable  to  the  right 
or  left. 

3  This  screw  adjusted  the  platform 

toward  or  away  from  the  viewer. 

All  of  the  above  functions  were  adjusted  in  a  specially  de¬ 
signed  fixture  using  precision  measuring  devices  that  mea¬ 
sured  to  0.  0001  inch. 

Figures  10c  and  lOd  show  the  fixture  being  used 
with  Discos  package  No.  5.  The  dial  indicator  (used  for 
the  measurement  of  the  vertical  motion  of  the  platform)  as 
well  as  the  height  gauge  (used  to  locate  the  package  counter¬ 
weight)  are  shown.  Figure  10  also  shows  the  micrometers 
used  for  the  lateral  positioning  of  the  package. 


Arrow 

1 

2 
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Fig.  7  PACKAGE  ON  BIFILAR  PENDULUM  WITH  ARROWS  SHOWING 
ADJUSTMENT  METHODS 
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Fig.  8  PACKAGE  ON  BIFILAR  PENDULUM  IN  ONE  OF  THE  ORTHOGONAL 
AXIS  POSITIONS 
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All  of  the  adjustments  were  made  after  the  center  of 
mass  was  measured  (see  Figs.  10a  and  10b).  The  beam 
used  for  the  center  of  mass  determination  was  designed  so 
that  the  moment  of  inertia  platform  was  registered  in  a 
known  location  on  the  beam.  This  beam  was  supported  by 
a  0.  375  inch  diameter  ball  in  the  center  of  the  weighing  pan 
on  a  precision  balance  and  two  0.  375  inch  diameter  balls  at 
the  ciner  end. 

The  advantage  of  measuring  the  center  of  mass  in 
this  manner  was  that  the  package  (once  installed  on  the 
inertia  platform)  was  not  removed  until  all  of  its  mass  prop¬ 
erties  were  determined. 

To  complete  the  system,  the  pendulum  was  sus¬ 
pended  in  a  bell  jar  (Fig.  11),  which  rested  on  a  collar. 

The  collar  had  two  ports  that  could  be  used  for  viewing, 
and  when  required  the  port  cover  could  be  disassembled 
for  system  adjustment  without  disturbing  the  bell  jar. 
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5.  ELECTRONICS  INSTRUMENTATION 


It  was  obvious  in  the  development  of  the  measuring 
technique  that  the  bifilar  pendulum  was  going  to  have  to  be 
carefully  instrumented  to  be  successful.  Figure  11  shows 
all  of  the  instruments  external  to  the  be’ l  jar,  and  Fig.  8 
is  a  view  that  shows  the  instrumentation  inside  the  bell  jar. 

A  shaft  was  attached  to  the  bottom  of  the  moment  of 
inertia  platform.  It  was  guided  by  a  nylon  bearing  approxi¬ 
mately  0.  010  inch  larger  than  its  O.  D.  To  start  the  pendu¬ 
lum,  a  timing  motor  was  actuated  that  rotated  a  cam- 
operated  wire  that  "locked"  onto  the  pendulum  shaft.  After 
the  pendulum  was  rotated  20°,  the  wire  was  abruptly  re¬ 
leased  from  the  inertia  pendulum  shaft,  allowing  the  pendu¬ 
lum  to  swing  freely. 

A  fan-shaped  disk  with  a  small  hole  in  it  was  located 
directly  above  the  pendulum  shaft.  The  "fan"  rotated  be¬ 
tween  a  light  and  an  N-P-N  planar  silicon  light  sensor. 

When  the  hole  passed  over  the  light  sensor,  the  light  illumi¬ 
nating  the  sensor  triggered  the  external  electronics. 

The  following  is  a  list  of  the  external  electronics 
and  their  use: 

Instrument  Function 

Oscilloscope  Tc  determine  if  the  light  was 

centered  over  the  hole  in  the 
"fan.  " 

Sanborn  recorder  To  check  for  pendulum  damp¬ 
ing. 

Printer  To  record  the  period  of  the 

pendulum. 

Counter  To  check  if  the  average  of  10 

periods  was  remaining  con¬ 
stant. 


Preceding  page  blank 
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6.  TEST  RESULTS 


On  23  November  1971  the  final  calibration  tests  and 
package  testing  began  on  the  bifilar  pendulum.  During  all 
of  the  tests  a  vacuum  of  30  inches  of  mercury  was  main¬ 
tained. 


To  account  for  any  measurement  variation  owing  to 
the  weight  of  the  package  being  measured,  a  calibration 
weight  was  made  for  each  package.  Also,  to  determine  if 
the  measured  moments  of  inertia  were  reasonable,  a  com¬ 
puter  program  was  written  that  calculated  the  moment  of 
inertia  of  homogeneous  packages  the  same  size  as  Discos 
(see  Appendix  A).  Figure  12  shows  the  program  printout. 

To  do  the  final  calculations  a  computer  program  was 
written.  This  program  calculated  the  test  mass  moments 
of  inertia  as  well  as  the  final  package  moment  of  inertia. 

In  this  way  the  program  was  continually  checked  for  accu¬ 
racy  by  a  known  mass.  Figure  13  is  a  sample  of  the  pro¬ 
gram  printout.  The  printout  notation  is  as  follows: 


I..P  =  Final  package  moment  of  inertia  and  cross 


U 


product. 


IL(ij)  =  Moment  of  inertia  of  skew  line, 
li  =  Principal  moment  of  inertia, 
i  =  ] ,  2,  3,  and 

IPRO(ij)  =  Final  moment  of  inertia  matrix  check. 
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Fig.  12  HOMOGENEOUS  PACKAGE  PRINTOUT 


Fig.  13  FINAL  PROGRAM  PRINTOUT 
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7.  CONCLUDING  REMARKS 


On  11  April  1972  Figs.  14a  and  14b  were  sent  to 
Stanford  University,  thereby  completing  the  package  mass 
property  documentation  for  the  Discos  electronic  packages. 

After  insertion  into  orbit.  Discos  operated  success¬ 
fully  for  1  year  until  it  was  commanded  off.  It  exceeded  its 
design  requirements  for  orbital  perturbations  of  10"H  g. 
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Appendix  A 

MOMENT  OF  INERTIA  EQUATIONS  FOR  BODIES  THE 
SHAPE  OF  THE  DISCOS  PACKAGES 


To  determine  if  the  moments  of  inertia  calculated  on 
the  bifilar  pendulum  were  reasonable,  the  equations  for  the 
moment  of  inertia  of  the  three  principal  axis  of  the  Discos 
package  were  determined  as  follows: 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  FHYSICS  LABORATORY 

Silvcr  Spuing,  Maryland 


The  radius  to  the  element  of  mass  is 


D  =  [z2  +  r2  cos2  6]  *  . 


(A-l) 


The  element  of  mass  is 


dm  =  Pr  dQdrdz  . 


(A-2) 


Therefore,  the  M  of  Ixx  is 


Ixx  =  p  J  JJ*  D  r  d0drdz  , 


H  R 


Ixx  =  p 


if } 

JK  •'r  -ot 
"2 


(A-3) 


(z2  +  r2  cos26!>  r2d0drdz  ,  (A-4) 


and,  finally. 


HpatR^  -  Rj]  Hp  sin2or  [R^  -  R*] 

IxX  = - 4 -  + - 8 - 


(A-5) 


H3ap  [R2  -R2] 
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And,  .finally, 

T  hpR4a  2hpR4  sin4a  hpR4sin4a  ,pR2h3a 
Iyy  4  is  2“tana  + - 12~ 

L 

f  pR2h3  sin2  a  pR2h3  sin2  a 
24  12  tana 

Izz 


The  element  of  mass  is 


dm  =  rdQdrdzp 


(i5 
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(A-13) 


The  above  equations  were  programmed,  and  the 
mass  properties  of  homogeneous  packages  of  the  same 
weight  as  the  actual  Discos  packages  were  determined. 
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